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ABSTRACT 

The  objective  of  this  study  is  to  examine  the  impact  of  single 
and  multi-component  surrogate  fuel  mixtures  on  the 
atomization  and  mixing  characteristics  of  non-reacting 
isothermal  diesel  engine  sprays.  An  Eulerian  modeling 
approach  was  adopted  to  simulate  both  the  internal  nozzle  flow 
dynamics  and  the  emerging  turbulent  spray  in  the  near  nozzle 
region  in  a  fully-coupled  manner.  The  Volume  of  Fluids  (VoF) 
methodology  was  utilized  to  treat  the  two -phase  flow  dynamics 
including  a  Homogenous  Relaxation  approach  to  account  for 
nozzle  cavitation  effects.  To  enable  accurate  simulations,  the 
nozzle  geometry  and  in-situ  multi-dimensional  needle  lift  and 
off-axis  motion  profiles  have  been  characterized  via  the  X-ray 
phase-contrast  technique  at  Argonne  National  Laboratory.  The 
flow  turbulence  is  treated  via  the  classical  k  —  e  Reynolds 
Average  Navier  Stoke  (RANS)  model  with  in-nozzle  and  near 
field  resolution  of  30  qm.  Several  multi-component  surrogate 
mixtures  were  implemented  using  linear  blending  rules  to 
examine  the  behavior  of  petroleum,  and  alternative  fuels 
including:  JP-8,  JP-5,  Hydro-treated  Renewable  Jet  (HRJ),  Iso- 
Paraffinic  Kerosene  (IPK)  with  comparison  to  single¬ 
component  n-dodecane  fuel  on  ECN  Spray  A  nozzle  spray 
dynamics.  The  results  were  validated  using  transient  rate-of- 
injection  measurements  from  the  Army  Research  Laboratory  at 
Spray  A  conditions  as  well  as  projected  density  fields  obtained 
from  the  line-of- sight  measurements  from  X-ray  radiography 
measurements  at  The  Advanced  Photon  Source  at  Argonne 
National  Laboratory.  The  conditions  correspond  to  injection 
pressure,  nominal  fuel  temperature,  and  ambient  density  of 
1500  bar,  363  K,  and  22.8  kg/w?,  respectively.  The  simulation 
results  provide  a  unique  high-fidelity  contribution  to  the  effects 
of  fuels  on  the  spray  mixing  dynamics.  The  results  can  lead  to 
improvements  in  fuel  mixture  distributions  enhancing 
performance  of  military  vehicles. 


Chol-Bum  M.  Kweon 

U.S.  Army  Research  Laboratory 
Aberdeen  Proving  Ground,  MD 


NOMENCLATURE 


M 

=  Projected  density^  ug/mm^ 

TIM 

-  Total  Integrated  Mass,  ug/mm 

do 

=  Nozzle  diameter,  um 

X 

=  Axial  distance  from  nozzle,  mm 

y 

=  Transverse  distance  from  nozzle,  mm 

Pf 

=  Density  of  liquid  fuel,  kg/m^ 

Pa 

-  Density  of  ambient  gas,  kg/m^ 

VOF 

=  Volume  of  Fluid  model 

a 

=  Volume  of  Fluid  Scalar 

ROI 

=  Rate  of  Injection,  mg/ms 

HRJ 

-  Hydro-treated  Renewable  Jet  Fuel 

IPK 

-  Iso-paraffinic  Kerosene  Jet  Fuel 

JP-8 

-  Jet  Propellant  8 

JP-5 

=  Jet  Propellant  5 

CFD 

=  Computational  Fluid  Dynamics 

RANS 

=  Reynolds  Average  Navier  Stokes 

ARL 

=  Army  Research  Laboratory 

AMR 

-  Adaptive  Mesh  Refinement 

MPI 

-  Message  Protocol  Interface 

PISO 

-  Pressure  Implicit  Split  Operator 

ANL 

-  Argonne  National  Laboratory 

ECN 

=  Engine  Combustion  Network 

EDM 

=  Electric  Discharge  Machining 

ASOI 

=  After  Start  of  Injection 
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INTRODUCTION 

An  important  factor  in  direct  injection  engines  is  characterizing 
the  spray  mixture  formation  process.  Improvements  in  the  air- 
fuel  mixing  can  lead  to  a  leaner  combustion  process  resulting  in 
enhanced  engine  efficiencies  and  improved  performance.  The 
spray  mixing  characteristics  can  be  strongly  influenced  by  the 
fuel’s  physical  properties  i.e.,  volatility,  viscosity,  and  surface 
tension.  These  properties  are  largely  dependent  on  the  fuel 
chemical  composition  and  the  hydrogen  to  carbon  ratio. 
Petroleum-based  fuels  are  chemically  complex  typically 
containing  thousands  of  components.  Large  variations  in 
petroleum  derived  fuel  composition  have  been  reported  arising 
from  individual  refinery  processes,  crude  oil  source,  and  also 
varying  with  season,  year  and  age  of  the  fuel.  This  myriad  of 
complex  factors  makes  it  difficult  to  control  the  consistency  of 
the  fuel  composition  for  research  purposes.  As  a  result, 
surrogate  fuels  are  desirable  for  experimental  and 
computational  tractability  and  reproducibility.  Surrogate  fuels 
have  typically  been  developed  to  represent  physical  or  chemical 
properties  and  are  particularly  useful  in  numerical  simulations 
since  it  simplifies  the  description  of  the  fuel  composition  while 
reducing  the  computational  cost  burden.  The  importance  of 
surrogate  fuels  in  research  (experiments  and  modeling)  has 
been  emphasized  through  the  launching  of  a  large  workshop 
named  Surrogate  Fuels  Working  Group  to  bring  forward  a 
consensus  on  the  development  of  databases  for  real 
transportation  fuels  (1).  As  a  result,  the  jet-fuel  (also  gasoline, 
and  diesel)  working  group  has  provided  a  palette  of  compounds 
from  which  to  construct  surrogate  fuel  mixtures;  several  of  the 
existing  surrogates  have  also  been  presented  and  compared  to 
experimental  data  (2-11). 

Previous  numerical  works  have  focused  on  developing  the 
chemical-kinetic  models  describing  select  aspects  of  kerosene 
spray  behavior.  Earlier  works  by  Schultz  et  al.  (2)  proposed  a 
12-component  surrogate  mixtures  which  has  led  to  further 
developments  of  mixtures  for  various  applications.  The  work 
by  Edwards  et  al.  (3)  have  reviewed  several  surrogate  fuels  and 
targeted  validations  with  various  applications.  It  was  noted  that 
fuel  injection,  vaporization,  and  mixing  should  use  a  multi- 
component  mixture  where  the  distillation  curve  is  a  critical 
parameter.  Bravo  et  al.  (4)  carried  out  high-fidelity  non-reacting 
spray  simulations  and  has  shown  the  suitability  of  several 
multi-component  kerosene  surrogates  at  diesel  engine 
conditions.  The  measurements  were  validated  with  the  Army 
Research  Laboratory  (ARL)  evaporating  spray  measurements 
of  Kurman  et  al.  (5)  demonstrating  good  agreement  for  several 
fuel  injector  configurations.  Under  reacting  conditions,  Violi  et 
al.  (6)  presented  a  surrogate  mixture  of  six  pure  hydrocarbon 
(Utah  surrogate)  and  found  that  it  correctly  simulated  the 
distillation  characteristics  of  JP-8.  Their  work  was  validated 
with  laminar  premixed  flames  burning  kerosene.  Similarly, 
Montgomery  et  al.  (7)  developed  a  reduced  four  component 
mixture  mechanism  with  reasonable  agreement  to  a  detailed  12 
component  surrogate  and  ignition  delay  measurements.  More 
recently,  Humer  et  al.  (8)  introduced  three  additional  surrogates 


composed  of  key  reference  fuels  denoted  as  surrogates  A,  B, 
and  C,  respectively.  The  counterflow  configuration  was 
employed  for  validation  with  experiments  where  the  predictions 
of  critical  conditions  of  extinction  and  auto-ignition  for 
surrogates  were  in  good  agreements  with  jet  fuel  data.  Note  that 
Honnet  et  al.  (9)  also  introduced  a  widely  used  kerosene 
surrogate  for  JP-8,  the  Aachen  surrogate.  This  two  component 
mixture  was  targeted  at  reproducing  critical  conditions  in 
addition  to  volume  soot  fraction  measured  in  laminar  non- 
premixed  flows.  This  also  provided  an  enhanced  chemical 
mechanism  for  the  combustion  process  and  was  validated  with 
a  wide  set  of  data  from  shock  tubes,  rapid  compression 
machines,  jet  stirred  reactors,  burner  stabilized  premixed  flames 
and  freely  propagating  premixed  flame.  Encouraging  results 
were  obtained  with  respect  to  the  non-premixed  combustion  of 
kerosene.  It  is  important  to  note  that  until  recently  surrogate 
mixtures  have  not  been  validated  in  detail  numerically  for  non¬ 
reacting  spray  applications  with  military  and  alternative  fuels  at 
diesel  like  engine  conditions. 

Recent  developments  in  X-ray  radiography  at  The  Advanced 
Photon  Source  at  Argonne  National  Laboratory  (ANL)  have 
provided  unique  measurements  of  the  spray  structure  in  the 
near  nozzle  field  (10-20  mm).  Kastengren  et  al.  (13)  presented 
measured  projected  mass  density  fields  in  diesel  type  jets  at 
various  operating  conditions  demonstrating  the  asymmetries  in 
fuel  mass  distributions.  The  asymmetries  were  reported  to  arise 
from  the  manufacturing  eccentricities  leading  to  imperfections 
in  the  injector  nozzle  geometry  (i.e.,  non-axial  nozzle  exit)  (14). 
Similar  studies  have  also  measured  the  internal  transient 
geometry  and  motion  of  the  injector  nozzle,  needle  valve  and 
hydraulics  (14).  These  findings  have  sparked  interest  in  spray 
modeling  work  which  couples  the  transient  needle  motion  with 
the  internal  flow  dynamics  and  ensuing  spray.  The  works  by 
Xue  et  al.  (15-16)  have  integrated  this  in  an  Eulerian  modeling 
approach  using  a  single-component  reference  fuel  and  reported 
on  the  effect  of  needle  wobble  on  spray  dynamics  at  engine  like 
conditions.  The  work  also  includes  validations  with  X-ray 
radiography  projected  density  fields  and  measured  mass -flow- 
rates  profiles  with  good  agreement  (16). 

In  this  work,  a  group  of  surrogate  fuels  have  been  selected  to 
investigate  the  mixing  properties  of  non-reacting  heavy  fuel 
sprays.  A  unique  Eulerian  3D  Computational  Eluid  Dynamics 
(CED)  spray  modeling  approach  is  adopted  to  model  the 
atomization  and  mixing  process.  A  Reynolds -Averaged  Navier 
Stokes  (RANS)  simulation  turbulence  modeling  approach  is 
prescribed  in  addition  to  the  single-fluid  mixture  model  based 
Volume  of  Eluids  (VOE)  approach  for  the  two-phase  flow 
analysis.  The  ARL  experimental  dataset  corresponding  to  fuel 
injector  analyzer  measurements,  is  utilized  for  validation 
purposes  on  the  nozzle  mass  flow  rates  i.e.,  Rate-of-Injection 
(ROI).  Also  note  that  the  surrogate  physical  properties  have 
been  implemented  for  the  numerical  studies.  Note  that  a 
principal  objective  of  the  current  work  is  to  establish  the 
suitability  of  proposed  surrogate  mixtures  that  have  the  similar 
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physical  and  chemical  properties  as  of  a  practical,  petroleum- 
derived  and  alternative  fuels,  in  particular,  the  mixing  jet 
behavior  at  diesel  like  conditions.  This  validation  is  necessary, 
as  it  will  serve  as  guidance  to  future  evaporating  non-reacting 
and  reacting  engine  spray  simulations. 

NUMERICAL  METHOD 

The  multi-dimensional  CONVERGE  CED  solver  (17)  has  been 
adopted  in  this  study  to  perform  the  detailed  Eulerian  spray 
simulations.  The  software  is  a  compressible  Navier  Stokes 
solver  which  is  based  on  a  first  order  predictor-corrector  time 
integration  scheme,  and  a  choice  of  second  or  higher  order 
finite  volume  schemes  for  spatial  discretization.  It  features  a 
non-staggered,  collocated,  computation  grid  framework 
utilizing  a  Rhie-Chow  interpolation  technique  to  avoid  spurious 
oscillations.  An  efficient  geometric  multi-grid  treatment  is  used 
to  solve  the  pressure  equation,  and  parallel  computing  is  based 
on  implementation  of  Message  Passing  Interface  (MPI) 
protocol.  It  provides  the  option  of  increasing  resolution  locally 
through  static  fixed-grid  embedding,  and  dynamically  through 
Adaptive  Mesh  Refinement  (AMR)  activated  through  user 
specified  criteria.  The  solver  also  provides  a  choice  between  a 
number  of  modeling  options  for  the  treatment  of  turbulence 
including  advanced  Large  Eddy  Simulation  (LES)  and 
Reynolds -Averaged  Navier-Stokes  (RANS)  models. 
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The  R  parameter  is  a  function  of  turbulence  model  constants  for 
the  RANS  model.  Their  values  in  this  study  are  =  1.60, 

Q2  =  1-92,  and  Q3  =  —1.0  . 


In  the  above  K,  f ,  /i^,  are  conductivity,  turbulent  diffusion, 

gas  temperature,  and  enthalpy  of  species  of  the  mixture. 


The  fuel  mass  fraction  description  follows. 


The  Eulerian  formulation  adopted  in  this  study  describes  the 
ambient  gas  and  liquid  fuel  as  a  single  fluid  mixture  using  the 
Eavre  Averaged  Navier  Stokes  Equations  while  adopting  a 
traditional  RANS  modeling  (18-19)  approach.  The 
compressible  system  of  transport  equations  for  mass,  and 
momentum  are  presented  here. 


dpUj 

dXj 


=  0 


(1) 


dpY  ^  dpu,Y  _  dpu/Y' 


(8) 


with  u[  and  T' representing  fluctuations  in  velocity  and  liquid 
fuel  mass  fraction,  and  pYevap  representing  an  evaporation 
source  term  which  is  negligible  for  the  present  non-evaporating 
conditions.  The  closure  for  the  liquid  mass  transport  term 
pu^'Y'  is  based  on  a  turbulent  gradient  flux  model  (18), 


dpUi  dpUiUj  dP  da^j  dXij 

dt  dXj  dXi  dXj  dXj 


(2) 


-'T^,  dY 


(9) 


where  cry  and  Ty  are  the  viscous  and  modeled  Reynolds  stress 
tensor.  The  viscous  stress  is  modeled  using  the  Boussinesq  eddy 
viscosity  method  with  5y  as  the  Kronecker  delta, 


Tj/  =  Pt 


dui  dUj  2  3ufc 
dxj  dXi 


A  Standard  k  —  e  RANS  model  is  used  to  as  follows. 


(4) 


Where  =  0.7.  The  averaged  density  of  the  immiscible  two 
phase  mixture  is  given  by. 


1 _Y  1-Y 
P  Pi  Pg 


(10) 


where  gas  density  is  prescribed  through  an  equation  of  state 
=  p/(Rgf)  while  assuming  the  liquid  phase  density  to  be 
constant.  It  is  important  to  note  that  multi -component  mixture 
properties  are  prescribed  through  weighted  average  property 
values  from  the  surrogate  blend  mass  fraction  distributions. 


With  the  turbulent  kinetic  energy,  dissipation  and  energy 
equations  as. 


The  VOE  Eulerian  method  used  for  the  two-phase  flow  analysis 
is  based  on  the  volume  fraction  scalar(a).  It  is  defined  as 
follows. 


3 


10  cell  is  filled  with  pure  liquid 
1  cell  is  filled  with  pure  gas 

Intermediate  values  between  0  <  a  <  1  represent  a  mixture  of 
liquid  and  gas.  In  the  current  study,  the  volume  fraction  is 
directly  extracted  from  the  fuel  mass  fraction  transport  equation 
(8)  and  the  bulk  value  of  mass-densities  as  follows, 


(1  - 

{l-Y)/^  +  Y/pi 


(11) 


For  time  advancement,  equations  (1)  and  (2)  are  solved 
numerically  using  a  classical  predictor-corrector  scheme  in 
which  the  velocity  field  is  first  integrated  using  the  Navier 
Stokes  equation  (2)  and  then  corrected  to  enforce  mass 
conservation  (1),  using  a  modified  pressure  value.  This  is 
achieved  using  the  modified  Pressure  Implicit  with  Splitting  of 
Operator  (PISO)  algorithm  first  introduced  by  Issa  (12). 


EXPERIMENTS  FOR  COMPARISON 

Two  datasets  of  measurements  were  utilized  in  this  study  for 
validation  purposes.  The  first  dataset  corresponds  to  Argonne 
National  Laboratory  (ANL)  X-ray  radiography  measurements 
of  spray  projected  density,  nozzle  geometry,  and  needle 
displacement  for  single-component  n-dodecane  fuel  (20).  Note 
the  projected  density  field  measurements  are  presented  as 
ensemble  averaged  fields  over  multiple  shots.  The  table  below 
shows  the  experimental  conditions: 

Table  1.  Conditions  for  non-evaporating  Spray  A  X-ray 
radiography  measurement  (21). _ 


Ambient  gas  temperature 

303  (K) 

Ambient  gas  pressure 

2.0  (MPa) 

Ambient  gas  density 

22.8  (kg/m^) 

Ambient  gas  N2 

100% 

Nozzle  K  factor 

1.5 

Nozzle  outlet  diameter 

0.090  mm 

Number  of  holes 

single-hole 

Fuel  injection  pressure 

150  (MPa) 

Fuel 

n-dodecane 

Fuel  temperature  at  nozzle 

343  K 

Injection  duration 

1 .5  (ms) 

Injected  mass 

3.5  (mg) 

The  second  dataset  was  conducted  at  the  ARL  using  a  common- 
rail  fuel  injector,  lAV  injection  analyzer,  for  nozzle  mass-flow- 
rates,  i.e.,  Rate-of-Injection,  measurements  of  JP-8  and  IPK 
fuels.  It  is  important  to  note  that  ROI  profiles  are  also  averaged 
over  100  shots  each  reducing  the  effects  of  outliers  in  the 
signal.  The  measurement  details  for  this  study  are  provided 
below. 


Table2.  Conditions  for  nozzle-flow  dataset  ROI  measurements 
with  various  fuels. 


Injector  type 

BOSCH  CR1N3 

Number  of  holes 

single-hole 

Fuels 

n-dodecane,  JP-8,  IPK 

Fuel  injection  pressure 

150  (MPa) 

Nozzle  K  factor  (nominal) 

1.5 

Nozzle  outlet  diameter 

0.090  mm 

Injection  duration  (nominal) 

1.5  (ms) 

COMPUTATIONAL  SETUP 

Three-dimensional  simulations  have  been  conducted  using  a 
single-mixture  VOF  model  that  couples  the  needle  motion  with 
the  ensuing  spray.  The  nozzle  tip  injector,  sac,  and  needle 
geometry  have  been  utilized  in  the  simulations  via  the 
characterization  obtained  from  the  Engine  Combustion 
Network  (ECN)  community  via  high-resolution  X-ray 
measurements.  Figure  1,  shows  the  characterization  of  needle 
valve  motion  obtained  from  (20). 


Time  (ms) 

Figure  1.  Spray- A  Injector  characterization.  Needle  lift  and  off- 
axis  measurements. 

The  transient  needle  lift  shown  along  the  y-axis  in  the  figure 
prescribes  most  of  the  injected  fuel  mass  in  the  ROI 
measurements.  The  off-axis  motion  results  in  needle  wobble  in 
the  transverse  and  span- wise  flow  directions.  The  needle 
profiles  are  prescribed  as  dynamic  boundary  conditions 
(immersed)  in  the  simulation  work. 

The  spray  chamber  region  in  the  simulations  (cylindrical)  has 
dimensions  of  50  mm  by  200  mm  in  diameter  and  length 
respectively.  Fixed  embedding  is  prescribed  in  the  near  nozzle 
region  with  a  base  grid  size  of  2  mm.  Figure  2a  shows  a  fuel 
mixture  distribution  for  IPK  at  t  =  0.51  ms,  highlighting  the 
nozzle  geometry.  Figure  2b  shows  a  close  inspection  of  the 
injector  nozzle  with  the  cell  distribution  and  fixed  embedding 
regions.  Figure  2b  also  shows  the  misalignment  of  the  orifice  in 
the  sac  arising  from  the  manufacturing  eccentricities. 
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Figure  2.  Coupled  spray  solution  and  mesh,  (a)  Fuel  mixture 
fraction  distribution,  and  (b)  zoomed  view  of  the  mesh  in 
nozzle  and  outflow  regions. 


From  previous  simulation  studies  it  was  concluded  that  a  grid- 
size  of  30  |im  (nozzle  diameter  to  cell  size  ratio  of  3)  would  be 
sufficient  to  capture  the  important  spray  features  (16).  This 
recommendation  is  based  on  a  trade-off  between  accuracy  and 
computational  cost.  Note  also,  each  simulation  was  run  on  32 
processors  with  peak  cell  counts  of  --850,000  and  total  wall- 
clock  time  of  128  hours  running  on  CRAY  XE6  Linux 
architectures.  The  resolution  is  embedded  inside  the  nozzle 
orifice  and  up  to  2  mm  axially  in  the  spray  region.  The 
surrogate  components  for  each  fuel  are  listed  in  Table  3  in 
molar  fraction. 

RESULTS  AND  DISCUSSIONS 

The  Eulerian  VOF  simulations  were  carried  out  and  inspected 
using  a  resolution  of  30  |im  in  the  nozzle  and  optically  dense 
spray  regions.  A  single  RANS  realization  of  the  mean  flow  field 
was  obtained  for  each  of  the  single  and  multi-component  fuels 
considered.  The  fuel  distributions  are  compared  using  the 
projected  mass  density  from  the  simulations.  Note  that 
experimentally  this  is  obtained  from  a  line-of-sight  integration 
along  the  x-Ray  beam.  To  enable  one-to-one  comparison  a 
similar  procedure  was  implemented  while  post-processing  the 
simulation  results. 

Table  3.  Surrogate  components  and  its  blend  density  at  300K 


TRANSIENT  RESULTS 

The  projected  mass  density  simulation  results  below  compare 
the  influence  of  density  variations  in  single  and  multi- 
component  surrogate  mixtures.  The  solution  corresponds  to  an 
early  integration  time  of  t  =  0.025  ms.  It  is  apparent  that  each 
surrogate  fuel  has  a  specific  and  characteristic  mixture 
distribution.  Note  the  peak  penetration  lengths  for  IPK  and  JP-8 
occur  at  8.8  mm  and  8.9  mm,  respectively  with  density  values 
of  M  =  25  lug/mm^.  Although  n-dodecane,  HRJ  and  JP-5  have 
similar  peak  penetration  values  -6.8  mm,  it  is  clear  that  the 
spray  structure  is  different  as  indicated  from  the  dispersion  and 
spray  tip  characteristics. 


X  (mm) 

Figure  3.  Contours  of  projected  mass  density  (M)  distributions 
(ug/mm^),  at  0.025  ms  AS 01  from  simulations  with  various  fuel 
mixtures:  (a)  n-dodecane  (b)  HRJ  (c)  IPK  (d)  JP-5  (e)  JP-8. 


Fuel 

Surrogate  component  (mole  fraction) 

Surrogate 

blend 

density 

(kg/m^ 

Hexadecane 

Iso- 

dodecane 

m-Xylene 

Butyl 

cyclohexane 

n-Propyl 

benzene 

Decane 

Trimethyl- 

Benzene 

JP5 

0.53 

0.15 

0.24 

0.08 

0.0 

0.0 

0.0 

790 

HRJ 

0.65 

0.33 

0.0 

0.0 

0.02 

0.0 

0.0 

762 

IPK 

0.21 

0.65 

0.0 

0.12 

0.02 

0.0 

0.0 

758 

JP8 

0.0 

0.0 

0.0 

0.0 

0.0 

0.83 

0.17 

756 
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Figure  4.  Transverse  profiles  of  projected  mass  density  (M)  at  4  axial  stations  {x  =  0.1,  2.0,  4.0,  6.0  mm)  in 
the  spray  downstream  direction  at  r  =  0.025  ms.  The  surrogate  mixtures  considered  correspond  to 
(a)  n-dodecane  (b)  HRJ  (c)  IPK  (d)  JP-5  and  (e)  JP-8. 
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Figure  4  presents  transverse  profile  of  the  projected  density  at 
several  locations  {v  =  0.1,  2.0,  4.0,  6.0  mm)  downstream  of  the 
nozzle.  It  is  apparent  that  each  transverse  profile  is  asymmetric 
arising  from  the  manufacturing  discrepancies  in  the  real  nozzle 
geometry.  The  projected  density  profiles  follow  a  typical  power 
law  M  -  MoX  "  accounting  for  the  decay  in  its  peak  magnitude 
with  axial  direction  (where  Mo  is  an  initial  value  and  ^  is  a 
power  constant).  The  jet  profiles  also  increase  in  width  with 
axial  direction  due  to  the  natural  increased  mixing,  i.e.,  air 
entrainment,  in  the  radial  direction.  Note  the  jet  width,  as 
indicated  by  the  zero  intercept  of  the  profiles  at  v  =6.0  mm,  has 
a  value  of  yjet 0.6  mm  for  IPK  and  JP-8  mixtures.  This  is  1.5 
times  larger  than  the  findings  for  n-dodecane,  HRJ,  and  JP-5 
with  similar  widths  at  yjet  0.44  mm  at  the  present  transient 
conditions.  The  results  suggest  that  differences  in  the  jet 
spreading  rate  arise  mainly  from  blend  density  variations  with 
lighter  fuels  typically  spreading  faster  axially  and  radically  than 
the  heavier  counterparts  (i.e.,  yjet(jp-8)  >  yjet(dodecane)). 


STEADY  RESULTS 

Figure  5  shows  steady  results  obtained  at  time  ^  =  0.51  ms  via 
contours  of  projected  mass  density.  The  simulations  are 
compared  against  the  X-Ray  n-dodecane  measurements. 


Figure  5, 
lug/mm^, 
mixtures 


4  6 

X  (mm) 

Contours  of  Projected  mass  density  (M)  distributions 
at  0.51  s  AS 01  from  simulations  with  various  fuel 
:  (a)  n-dodecane  (b)  HRJ  (c)  IPK  (d)  JP-5  (e)  JP-8. 


At  steady  conditions.  Figure  5  shows  insignificant  variations 
with  fuel  mixtures.  Near  the  axial  region  upstream  x  =  0-3  mm 
the  density  contours  of  M  show  undulations  induced  from  the 
off-axis  needle  motion.  Next,  the  full  spray  projected  density 
field  will  be  utilized  to  investigate  the  Total  Integrated  Mass 
(TIM)  variations  obtained  through  transverse  integration  of 
projected  density  fields.  The  projected  velocity  is  calculated  as 
the  inverse  of  TIM,  from  continuity  arguments,  and  normalized 
by  the  peak  injection  velocity  at  the  nozzle  exit. 


Axial  Direction  (mm) 


Figure  6.  Comparison  with  X-ray  measurements  at  steady 
conditions  t  =  0.51  ms  for  each  fuel  considered,  (a)  projected 
density  in  the  axial  direction,  [ug/mm^  (b)  Transverse  Integrated 
Mass,  TIM,  jug/mm  (c)  Relative  Velocity. 
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Figure  7.  Transverse  profiles  of  projected  mass  density  (M)  at  4  axial  locations  {x  =  0.1,  2.0,  4.0,  6.0  mm)  in 
the  spray  downstream  direction  at  r  =  0.51  ms.  The  surrogate  mixtures  considered  correspond  to  (a)  n-dodecane 
(b)  HRJ  (c)  IPK  (d)  JP-5  and  (e)  JP-8. 
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Figure  6a  shows  the  axial  development  of  centerline  projected 
density  including  the  comparison  between  n-dodecane  and  the 
measurements.  The  model  is  able  to  capture  the  correct  trend  as 
in  the  measurements;  however,  it  is  pointed  out  that  an 
optimization  of  the  parameters  (i.e.,  Schmidt  number  and 
turbulence  model  constants)  and  higher  resolution  may  be 
necessary  to  improve  the  results.  Figure  6a  profiles  show  that 
lighter  multi-component  mixture  such  as  JP-8  feature  increases 
in  projected  density  values  than  the  single-component  surrogate 
(n-dodecane)  due  to  their  faster  spreading  rates.  The 
comparisons  between  Figure  6b -6c  with  measurements  also 
show  good  agreements.  The  variations  with  different  fuels 
remain  subtle;  however,  it  is  evident  that  TIM  profiles  are 
larger  for  single-component  heavier  fuels  decreasing  with 
lighter  multi-component  mixtures  such  as  JP-8.  Similar 
findings  are  shown  with  the  relative  velocity.  Figure  7  depicts 
the  variations  of  transverse  projected  density  profile  with 
downstream  distance  {x  =  0.1,  2.0,  4.0,  6.0  mm}  of  the  spray  at 
r  =  0.51  ms.  The  nominal  jet  width  across  all  fuel  blends  is  jjei  - 
0.82  mm,  in  comparison  transient  profiles  report  a  peak  increase 
of  1.34  in  jet  width.  Note  that  the  variations  across  all  presented 
fuels  will  be  amplified  at  higher  ambient  temperatures  due  to 
differences  in  volatility  and  vapor  pressures  of  each  fuel. 


trend  in  peak-injected  mass  following  the  same  behavior  as  the 
CFD  counterpart.  Figure  8  (b)  is  a  scaled  up  figure  of  the  ROI 
simulations  showing  the  impact  of  physical  property  variations, 
in  particular  the  blend  fuel  density.  The  trend  from  larger  to 
smaller  ROI  values  becomes  apparent  in  this  figure  as  the 
following:  JP8>  IPK>  JP5>  HRJ>  dodecane  following  the 
surrogate  fuel  blend  prescribed  density  distributions  and  in 
agreements  with  experiments. 
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NOZZLE  MASS  FLOW  RATES  (ROh 

Measured  rate-of-injection  (ROI)  profiles  were  used  as  bulk 
validation  parameters  for  the  nozzle  flow  simulation.  As 
discussed  earlier,  the  fuel  injection  system  consists  of  a  high- 
pressure  fuel  bench  connected  to  a  common  rail  system.  The 
high-pressure  fuel  bench  contains  a  fuel  pump  operated  by  an 
electric  motor.  For  the  experiments  presented,  fuel  rail  pressure 
was  held  constant  at  1500  bar  following  Spray- A  specifications. 
Table  2.  The  motor  speed  is  maintained  by  a  variable  frequency 
drive.  Total  energizing  times  for  the  presented  experiments 
were  770  \xs  leading  to  an  injection  duration  of  1.5  ms.  For 
the  single-hole  ART  nozzle,  all  stock  6-holes  were  welded  and 
an  axial  orifice  was  created  by  an  electrical  discharging 
machine  (EDM)  machining.  Each  fuel  injector  was  then 
mapped  with  an  lAV  fuel  injection  analyzer  to  determine  the 
injected  fuel  mass  and  ROI  parameters  of  interest.  Note  that  the 
analyzer  operates  on  the  principle  that  the  fuel  mass  is  related 
to  the  speed  of  sound  in  the  fluid,  the  cross  sectional  area  of  the 
tube  and  the  wave  dynamic  pressure  as  a  function  of  time.  ROI 
signals  are  ensemble-averaged  over  100  realizations  to  reduce 
the  effects  of  instantaneous  fluctuations. 

Figure  8  (a)  shows  similar  behavior  across  all  measured  mass 
flow  rates.  However,  the  initial  ramping  of  measured  ROI 
differs  between  single  and  multi-component  fuels.  These 
variations  may  be  due  to  the  different  injector  types  used  in  the 
experiments  as  the  n-dodecane  measurements  correspond  to 
CRIN2  injector  (note  SNL  measurement  technique  is  different 
from  the  ARL  injection  analyzer.  SNL’s  ROI  data  overestimates 
ROI).  An  indication  of  this  is  the  consistency  between  IPK  and 
JP-8  dataset  in  the  ramping  region.  For  all  cases,  the  peak  ROI 
values  at  steady  conditions  (t  =  0.51  ms)  shows  the  increasing 


Figure  8.  Rate-of-Injection  profiles  at  r  =  0.51  ms  ASOI  from 
simulations  with  various  fuel  mixtures,  (a)  Comparison  of  with 
experiments  (b)  Highlights  of  variations  in  CFD  ROI 
distributions. 

At  initial  time  the  simulations  have  a  finite  needle  lift  of  20  um 
and  are  started  also  with  liquid  filled  sac  and  nozzle  hole.  This 
is  different  from  the  measurements  that  may  have  gas  in  sac 
and  nozzle  flow.  Also  observe  that  ROI  simulation  signals  show 
evidence  of  oscillations  across  all  fuels.  This  can  be  a  result  of 
presenting  single-realization  simulation,  as  some  of  the 
fluctuations  will  be  damped  through  averaging  of  additional 
events. 

CONCLUSIONS 

Three-dimensional  numerical  simulations  of  engine  sprays  have 
been  conducted  using  a  compressible  VOF  single  mixture 
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formulation  to  treat  the  multiphase  flow.  The  simulation 
geometry  includes  transient  needle  displacements,  sac  volume, 
and  the  nozzle  tapered  channel.  The  asymmetries  are  evident  in 
the  geometry  as  was  demonstrated  in  Figure  2.  The  fuels 
considered  in  this  study  correspond  to  practical  petroleum  and 
alternative  fuels  including  n-dodecane,  HRJ,  IPK,  JP-5  and  JP- 
8.  Single  and  multi-component  surrogate  mixtures  have  been 
selected  for  their  representation  as  presented  earlier.  The 
following  observations  are  summarized: 

•  RANS  simulation  results  are  adequate  in  capturing  the 
spray  injector  transients  in  jet  width  dispersion  and 
fuel  mass  density  distributions.  The  simulations 
compare  well  with  the  nozzle  and  ensuing  spray 
presented  measurements. 

•  Transient  results  of  projected  density  indicate  that 
spreading  rate  is  highly  depended  on  fuel  blend  mass 
density  distributions.  Jet  widths  for  IPK  and  JP-8  iyjet'- 
0.61  mm)  mixture  were  found  to  be  1.5  times  larger 
than  n-dodecane,  HRJ,  and  JP-5  iyjet-  0.44  mm). 

•  The  differences  in  steady  state  results  of  projected 
density  are  less  pronounced  than  their  transient 
counterpart  for  different  fuels.  The  nominal  jet  width 
across  fuel  blends  (yjet  -  0.82  mm)  was  reported  to  be 
1.3  times  larger  than  peak  transient  profiles  {jjet  -  0.61 
mm).  At  steady  conditions,  variations  in  projected 
density,  TIM  and  relative  velocity  consistently  show 
the  influence  of  higher  blend  density  on  solution. 

•  ROI  simulations  capture  the  correct  physical  behavior 
when  compared  to  experiments.  The  simulations  show 
higher  injection  rates  for  lighter  fuel  blends  with  good 
agreements.  ROI  signals  also  show  natural  oscillations 
arising  from  traveling  pressure  waves  inside  the  nozzle 
flow. 

The  results  present  successful  simulations  of  internal  nozzle 
flow  in  a  single-hole  injector  with  attention  to  surrogate 
mixtures.  Future  works  will  include  extending  the  present  set  of 
results  to  higher  ambient  temperature  and  pressures  (typical 
engine  operating  ambient  conditions  are  850-950  K  and  40-80 
bar)  while  keeping  the  fuel-to-air  density  ratio  constant.  It  is 
expected  that  evaporation  effects  will  have  a  strong  impact  on 
projected  density,  dispersion,  and  jet  width  characteristics. 
Since  X-ray  diagnostics  are  not  yet  performed  under  high- 
temperature  conditions  and  optical  diagnostics  cannot  provide 
the  above  mentioned  details,  CFD  simulations  may  be  the  key 
tool  in  further  understanding  differences  in  mixing  and 
combustion  characteristics  of  single  and  multi -component  fuels. 
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